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Abstract 

A new method of estimating higher order perturbative coefficients is dis- 
cussed. It exploits the rapid, asymptotic growth of perturbative coefficients 
and the information on the singularities in the complex Borel plane. A com- 
parison with other methods is made in several Quantum Chromodynamics 
(QCD) expansions. 
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The ordinary perturbative expansions in weak coupling constant in quantum field theories 
are generally asymptotic expansions, with their perturbative coefficients growing factorially 
at large orders [1]. In practice, in many cases the asymptotic growth sets in quite early in 
perturbation, and the rapid growth of the coefficients becomes apparent already at first few 
orders. 

The asymptotic divergence of the perturbative coefficients implies singularities in the 
Borel plane. There are three kinds of known singularities, all on the real axis: those 
instanton-induced, ultraviolet renormalons, and infrared renormalons. 

In this paper, we show that this rapid growth of the perturbative coefficients and the 
information on the singularities in the Borel plane can be turned into a useful tool that 
allows us to estimate the unknown (N + l)th order coefficient using the known coefficients 
up to order N. The method we are going to present can be most easily understood by 
working out an explicit example, the double-well potential in quantum mechanics. The 
tunneling between the two potential wells splits the degenerate perturbative ground states 
into a parity even, and an odd states, and the average energy E{a) of the energies of these 
two states has the perturbative expansion of the form 



E{a) = 



n=l 



(1) 



where a denotes the canonical coupling of the model [2]. The Borel transform E{b) of E{a), 
which has the perturbative expansion 



~ a 



1 + 



n=l 



(2) 



is known to have multi-instanton-anti-instanton caused singularities at & = 2nSo, {n — 
1, 2, 3, • • •), where yS'o = 1/6 is the one-instanton action (in units of 1/a). 

The nature of the singularities can in principle be determined by doing perturbation in 
the background of multi-instanton-anti-instanton configurations. The closest singularity to 
the origin at 6 = 1/3, which determines the leading large order behavior of the expansion 
(1), can be shown to have the following form [3] 



9 

m 



7r(l - 3by 



53 

1 - — (1 - 36) + 0[(1 - 36)2 ln(l - 36)] 

18 



-I- Analytic part , (3) 



which comes from the instanton-anti-instanton contributions to E{a). The "Analytic part" 
denotes terms that are analytic around the singularity. 

We now consider the function R{b) = (1 — 2>bYE{b) as introduced in [4], to control the 
divergence at the singularity. R{b) is bounded at 6 = 1/3, and has a very soft singularity, a 
suppressed logarithmic cut. Thus the power expansion of R{b) around the origin is expected 
to better behave than that of E{b). Ignoring the residual logarithmic cut, we would expect 
the convergence radius of the former is bounded by the second singularity at 6 = 2/3, thus 
effectively to become twice that of the latter. To make the expansion of R{b) better, we can 
take a further step of conformally mapping the singularities except the first one as far away 
as possible from the origin. This way one can hope to have a smoother R in the new plane, 
and a better behavior of the expansion. 
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Double well-potential 





ai 


02 


as 


04 


as 


Exact 


4.5 


44.50 


626.625 


1.103 10"^ 


2.288 10^ 


Est. 


6.0 


33.75 


600.047 


1.070 10^ 


2.266 10^ 



TABLE I. The exact and the estimated coefficients of the vacuum energy of the double-well potential. 

The conformal mappings we consider in this paper are such that all the singularities in 
the Borel plane except the first one are mapped to the unit circle and the first one to a point 
within the circle [5]. Such a mapping in this case is given by 

1 - Jl - 36/2 

w = \ (4) 

1 + -^1 - 36/2 

which maps the first singularity to Wq = (\/2 — l)/(v^ + 1) = 0.171 and aU others to the 
unit circle. Without the residual cut-singularity, R would be analytic on the unit disc in the 
w-plane. 

Now the power expansion of R{b{w)) up to 0{w^) is 



w 



n=l 



= 1 + (-16 + 2.667ai)w + (-120 + 3.55602)^^ + (-1896.44 + 3.161a3)w^ + 

(-22544 + 2.10704)^^ + (-254667.7 + 1.124a5)w^ + • • • (5) 
= 1 - 4w + 38.22^i;^ + 8Aw^ + 699.1w^ + 2538.1w^ + • • • , (6) 

where the coefficient at a given order in (5) was calculated using the exact values of a„, which 
are given in Ref . [2] , up to that order less one, and the coefficients in (6) were obtained with 
the exact values of a„. What is interesting about this expansion is that the a„-independent 
constant term in the coefficient at a given order in (5) is much larger than the exact value 
of the corresponding r„ in (6). This then implies that a good approximate estimate of a„ 
can be obtained by simply putting 

rn = (7) 

in (5). The estimated values for a„ from this prescription are given in Table I, which shows 
improving accuracy as the order grows, and the accuracy becomes better than 1% at n = 5. 

Obviously, the success of this method relies on the exact value of being much smaller 
than the constant term in the corresponding coefficient in expansion (5). We may understand 
this feature in the following way. We note first the coefficient r„ is a linear combination of 
aj's {i < n), which are large numbers. For the sake of argument, let us ignore for the moment 
the soft cut-singularity of i? at if = Wq. Then R is analytic on the unit disk in w-plane, 
and so we expect the growth of r„ is fundamentally limited. To yield a small number out of 
large numbers, this then suggests a delicate cancellation occur among the large numbers in 
the expression for r„. However, when r„'s are written as in expansion (5), the cancellations 
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are incomplete, yielding the large constant terms. This account shows that the essential 
ingredients for the success of our prescription are: (i) rapid growth of the perturbative 
coefficients, and (ii) information on the nature and locations of the Borel plane singularities. 
One may notice the exact values of r„'s in (6) still grow quite rapidly. This may be due to 
the residual logarithmic cut, and the infinite number of singularities on the unit circle. 

With this example in mind, we now give the general prescription for estimating uncalcu- 
lated higher order coefficients. Let us suppose an amplitude A[a) and its Borel transform 
in a given theory have the perturbative expansion of the form 

oo 

A{a) ^a+J2 dna''^^ (8) 

n=l 

with a denoting the canonical coupling of the theory, and 

oo Ln 

A{h) = 1 + ^ a„- . (9) 

n=l 

The overall normalization of the expansion (8) was chosen such that the coefficient of the 
leading term is unit. We assume the Borel transform A{h) has singularities on the real axis 
on the 6-plane, and the first two singularities on the positive real axis are at 6 = and at 
b = bi, respectively, and the first singularity on the negative real axis is at 6 = —b-Q. We 
further assume the nature of the first singularity on the positive real axis is known, and is 
of the form 

^(^) = + 0(1 - b/bo)] + Analytic part , (10) 

with C a real constant and the 'Analytic part' denoting terms analytic about b = bo. We 
now introduce a conformal mapping 



^1 + b/b^o - ^1 - b/b, 
w — — ( 11 ) 

^1 + + - V^*! ' 

which maps the first singularity on the positive axis to a point within the unit circle and all 
others to the circle, and the function 

R{b)^{l-b/bof+-A{b). (12) 

We then expand R{b{w)) in power series in w-plane: 

oo 

n=l 
oo 



n=l 



Note that g„ depend linearly on the leading coefficients to a„ only. The estimated nth 
order coefficient is then given by 
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Bjorken polarized sum rule 





J ext. 
Ul 


J est. 


«2 




U2 


«2 


//„cst. 


"3 




1X4 


1 


4.25 


3.86 


34.01 


27.36 


28.41 


27.25 


302 


290 


3750 


2716 


2 


3.92 


3.51 


26.94 


23.17 


24.09 


23.11 


223 


203 


2515 


1696 


3 


3.58 


3.14 


20.22 


19.21 


20.01 


19.22 


156 


130 


1565 


933 


4 


3.25 


2.73 


13.85 


15.46 


16.16 


15.57 


101 


68 


867 


396 


5 


2.92 


2.29 


7.84 


11.90 


12.59 


12.19 


60 


18 


388 


56 


6 


2.58 


1.79 


2.19 


8.48 


9.29 


9.08 


31 


-22 


95 


-115 


Adler function 


Nf 




















d4^^' 


1 


1.87 


3.45 


14.11 


8.18 


8.71 


7.54 


66 


75 


693 


550 


2 


1.76 


3.19 


10.16 


7.19 


7.55 


6.57 


43 


50 


391 


316 


3 


1.64 


2.90 


6.37 


6.26 


6.40 


5.61 


24.4 


27.5 


165 


151 


4 


1.53 


2.58 


2.76 


5.37 


5.27 


4.68 


9.9 


8.4 


8.04 


49 


5 


1.41 


2.22 


-0.69 


4.50 


4.16 


3.77 


0.35 


-7.7 


-87 


2.47 


6 


1.30 


1.80 


-3.96 


3.59 


3.08 


2.88 


-3.48 


-21 


-129 


4.18 



TABLE II. The Bjorken polarized sum rule and the Adler function: Estimated values along with the 
exact values and the PMS, ECH values when available. 



an = . (14) 

Qn 

The fractional error of this estimate is 

(15) 

and this scheme would work as long as r„ is much smaller than p„. Of course, it is impossible 
to know the error a priori, and the reliability of the estimated coefficients should be decided 
by other circumstantial informations, for instance, as in the following examples, by the 
pattern of the known terms in expansion (13), or by comparison with the estimated values 
from using different methods. Since the expansion (13) is in general not convergent on the 
whole unit disc due to the residual soft singularity, and because of the singularities on the 
boundary of the disc, it is possible in principle for r„ at a certain order to jump to a large 
value. In such a case the scheme could fail, but that does not seem to occur, at least, at the 
orders considered in the following QCD expansions. 

We now apply our scheme to the various QCD expansions, estimating the uncalculated 
next higher order coefficients, and compare the obtained values with the estimates from 
other methods. Of course, the coefficients are renormalization scheme and scale dependent, 
but for simplicity, in all of the following examples our consideration will be in the MS scheme 
and at the renormalization scale fi = Q, where Q denotes the energy scale of the problem in 
consideration. The canonical coupling in these examples is a{Q) = as{Q)/n, where as{Q) 
is the strong coupling constant. 

The Bjorken sum rule: Our first example is the QCD correction A for the Bjorken 
polarized sum rule, which has the perturbative expansion 
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oo 

A(a) = a + X] '^nd 

n=l 



(16) 



Incidently, this correction coincides with that of Gross-Lewellyn Smith (GLS) sum rule at 
next-leading order, and differs only by the small 'light-by- light' contribution at next-next- 
leading order. The first two coefficients di,d2 are known [6]. The locations of the first 
singularities of the Borel transform A (6) and the parameter u are given as [7] 

^0 = -^, 6i = |-, 6-0 = 7^, u=if3Jf3o-7o)/f3o, (17) 
Po Po Po 

where /3o, Pi are the first two coefficients of the QCD /3-function: 

(3o = (11 - 2/37V^)/4, /3i = (102 - 387V^/3)/16 , (18) 

and 7o is the one-loop anomalous dimension of the twist-four operator appearing in operator 
product expansion of the sum rule [8] : 

70 = (iVe-l/iVe)/3 = 8/9. (19) 

Here Nf, Nc denote the number of quark flavors and colors, respectively. These are all we 
need to estimate the unknown next higher order coefficients. 

The result is given in Table II, alongside of the estimates from the principle of min- 
imal sensitivity (PMS) method (and the effective charge (ECH) method) by Kataev and 
Starshenko [9]. The latter is based on the optimization of the perturbative amplitude over 
the parameter space of the renormalization scheme and scale dependence. It estimates the 
uncalculated higher order coefficients by reexpanding the optimized amplitude in terms of 
the coupling in the particular scheme and scale of interest. The recent, exact, partial cal- 
culation of the next-next-next-leading order correction [10] indicates this method, among 
others [11], works particularly well . 

The estimates from the two methods are in good qualitative agreement, with some ex- 
ceptions at d^. It is remarkable that these two completely different approaches offer such 
close estimates. As in PMS case, our method works best at Nf = 3, and the estimates 
become worse as Nf increases. Although the two approaches are in an overall agreement, 
the differences at ^3 estimates are large enough to be phenomenologically significant. For 
instance, the new estimate ^3 ^ 156 at Nf = 3, which was used in the recent GLS sum rule 
analysis [12], results in a less renormalization scale dependent amplitude than with the PMS 
estimate ^3 ~ 130. It is instructive to see this estimate more closely. The expansion of R, 
following (13), at Nf = 3, reads 

n=l n=l 

= 1 + (-3.72 + 1.185rfi)w + (-13.488 + 0.702^2)^^ + (-43.257 + 0.2774)w^ + • • • (20) 
= 1 + 0.527W + 0.709^2 + ■■■ (21) 

where the last two terms in (21) are exact. Notice that the first two p„'s are large compared to 
the corresponding r„'s, satisfying one of the absolutely necessary conditions for our scheme. 
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Quark mass 











j,_l^lar. /3o 






^2 res. 




est . 








1 


9.30 


7.33 


6.29 


12.1 


123.3 


118.8 


120.2 


117.6 


2225 


2248 


1894 


51428 


2 


8.52 


6.82 


5.83 


11.3 


104.8 


101.3 


102.3 


102.9 


1762 


1776 


1551 


37986 


3 


7.74 


6.28 


5.34 


10.5 


87.2 


84.9 


85.7 


89.2 


1357 


1364 


1251 


27098 


4 


6.96 


5.71 


4.83 


9.8 


70.7 


69.8 


70.6 


76.5 


1008 


1009 


993 


18476 


5 


6.18 


5.09 


4.28 


9.0 


55.1 


55.8 


57.2 


64.7 


712 


711 


774 


11848 


6 


5.40 


4.43 


3.76 


8.2 


40.5 


43.1 


45.5 


54.0 


469 


468 


589 


6959 



TABLE III. Estimated values for the quark mass expansion along with the exact values and the values 
from the residue, large (3o methods. 



and the pattern of the first two known coefficients suggests that pa is also hkely to be much 
larger than r^. If we believe that the first two r„'s have any indication on r^, it seems to be 
safe to assume \rs\ < 2.0. This then leads to 149 < da < 163. 

The estimates of the next higher order coefficients are obtained using the ^3 esti- 
mates. Compared to the PMS estimates there is some significant difference at this order. 
Since the obtained values depend on the estimated o?3, one should note that there is another 
uncertainty arising from the error in ^3 estimate. 

The Adler function: The Adler function D{a) for the vector current-current correlation 
function of different quark flavors has the perturbative expansion 

00 

D{a) ^a + J2 dna""^^ (22) 

n=l 

where the first two coefficients are known [13]. Several physical observables can be related 
to this function through the dispersion relations, and the perturbative coefficients of the 
former can be obtained once the corresponding coefficients in the Adler function are known. 
The locations of the first renormalon singularities and v in this case are given by [14] 

&o = |-, = &-o = -^, ^ = 2/3l//3o^ (23) 

Po Po Po 

and the resulting estimates for the first coefficients are given in Table II, along with the 
PMS, ECH values. Again, our method works best at iVj = 3 for the first two coefficients, 
and in this case the predicted 0^3 = 24.4 is slightly smaller than the PMS value 0^3 = 27.5. 

The quark mass: This last example is concerned with the on-shell quark mass. The 
on-shell quark mass and the MS mass are related by 



m-y^ 3 



n=l 



(24) 



and the first two coefficients are known [15]. The relevant parameters in this case are [16] 

&o = ^, &i = 4-, ^-0 = 7^, i^ = Am^ (25) 
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and the estimated first coefficients are in Table III, along with the values from the residue 
based method of Pineda [17] (also, for comparison, those of large Pq approximation [18]), 
which is known to work well in this case. The latter method relies on the rapid convergence 
of the perturbative calculation of the renormalon residue, and the expansion of the Borel 
transform about the first renormalon [4,5]. The agreement between these two approaches at 
rs estimates is remarkable. 

To conclude the paper, we presented a new method of estimating higher order unknown 
coefficients. It is based on the two generic features of the asymptotic expansions, the rapid 
growth of the coefficients and the Borel plane singularities. It can provide independent 
estimates for the yet unknown coefficients, which may prove useful in physical analysis. 
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Postscript: After submission of our paper we were informed that in the published 
version [19] of hep-ph/9811367 Borel transform with a conformal mapping was tried on a 
model function to estimate higher order coefficients. Unlike our method, however, it did not 
remove the first Borel singularity, and because of that it reached to a conclusion substantially 
different from ours, that the method can work only for large orders (say, N 7) and is not 
applicable for low orders (N~ 3) relevant for present QCD perturbations. 
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